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TOWARD AUTOMATED ANALYSIS OF PARTICLE HOLOGRAMS 


H.J. Caulfield, Aerodyne Research, Inc., 45 Manning Road, Billerica MA 01821 

ABSTRACT 

A preliminary study of approaches for extracting and analyzing data from 
particle holograms concludes that 

o For "thin” spherical particles out-of-focus methods are optimum, 
o For "thin" nonspherlcai particles out-of-focus methods are useful 
but must be supplemented by In-focus methods, and 
o A complex method of projection and back projection can remove the 
unwanted out-of-focus data for "deep" particles. 


I . INTRODUCTION 

Particle holograms are capable of Instantaneous recording of a large 
three-dimensional (3D) field particles. Given a particle hologram, we might 
have two tasks: 

1. Find all of the particles and 

2. Characterize them. 

The verbs "find" and "characterize" aust be defined n»re carefully for each 
particular situation. In most cases this analysis is carried out using a 
human observer to 

1. Locate the plane of best focus for each particle, 

2. Record the in focus image, 

3. Characterize that Image. 
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There are two major difficulties with this human-opera tor- based analysis: 
one obvious and one slightly more subtle. The obvious difficulty Is that 
human-based analysis Is notoriously slow, nonre pea table, and (hence) costly 
and inaccurate. It is this analysis bottleneck which keeps particle 
holography equipment on the shelf while far less powerful technologies perform 
a few of the tasks holography could perform. The less obvious difficulty is 
that human-based analysis may be Inherently less accurate than analysis based 
on a computer oriented algorithm. 

This paper explores past and current approaches to computer analysis of 
fields produced by particle holograms based on "multiplane** algorithms. 
Multiplane algorithms are based on data obtained in several depth planes 
(rather than simply the "focal** plane). 


II. PHILOSOPHICAL BACKGROUND 

A human observer looking at a point of laser light will see that point 
clearly. A recording surface moving through the field produced by such a 
point will produce a point image in focus and a diffraction pattern out of 
focus. This assymetry between the human (who sees only what is in focus) and 
the camera (which sees the out of focus data as well) leads to different data 
analysis schemes. In principle, all of the information in the particle 
wavefront is contained in every plane (including, of course, the hologram 
plane). Nonholographic cameras do not record the whole information, however, 
so observations in different planes give different information. 



To these fundamental observations we must add some practical ones. 
Particle hologram wavefronts are noisy and complicated. Other particles are 
present. Refractive artifacts occur. These complications mean that 
otherwise-equivalent observations are not really equivalent. More 
observations can mean reduced Influence of these "noise" effects. 


III. ANALYSIS OF SPHERICAL PARTICLES 


Spherical psrticlvS 
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known spheres only two questions can be asked: 
1. Where are they? 


and 


2- What are their diameters? 

>^£'/Che old (human-oriented) approach Is to move to the focal plane and 
there measure the lateral (x-y) position of the particle's center and measure 
its diameter. The depth (z) dimension Is that of the camera In the focal 
position. 

In the newer (computer oriented) approach, all of these quantities are 
measured In out of focus planes. We will review this past work briefly here. 

The first work in this field was by Vlkram and Billet. 1 They showed that 
diameter determination was far more accurate out of focus than in focus simply 
because the pattern la. In effect, magnified. Two derivative observations can 
be made immediately: 
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1. Localized film or detector noise is less of a problem out of focus 
than in focus and 

2. The centroid of the sphere (x— y location) can be located more 
accurately out of focus than in focus. 

The second work in this area was by Stanton, et.al. 2 They showed that 
1. The depth of the particle could be found more accurately by out of 
focus measurements than by in focus measurements and 
2* This obviates the need for searching all depth planes. 

The explanation of the increased accuracy in depth (z) location is a 
fairly universal one worthy of a little more explanation. If we plot "spot" 
size, s, versus z (z = 0 in focus), we obtain a curve with a minimum at z ■ 0. 
We can seek to find z ■ 0 by the z finding minimum s. Or we can measure s at 
two z * 0 positions on the same side of z « 0 and extrapolate to s ■ 0. The z 
sensitivity 

a “ 

is a maximum away from z ■ 0 and minimum at z ■ 0. Indeed maximum a occurs at 
maximum numerical aperture and way from z » 0. 

IV. SHALLOW OBJECTS 

A shallow object is one with unresolvable depth information. All that 
counts is its two-dimensional cross section normal in x-y. The questions 
which can be asked include 




1. What Is the x-y-z centroid location? 

2. What is the shape? (Not a well defined question) and 

3. What is the orientation? 

Whether these questions can be answered out of focus depends on what we mean 
by ’’shape" and "orientation." Suppose we mean by "shape" the best fit 
enclosing rectangle and by "orientation," the direction of the rectangle's 
predominant direction. Then clearly these are accessible out of focus. On 
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e.g. a particle shaped like a "6" and a particle shaped like a "9". 

When details of particle shape are of interest, we must go to a focal 
plane. This does not mean, however, we should Ignore the out of focus 
information. Extrapolation from out of focus data is still probably the best 
way to obtain focus. Furthermore, the out of focus data can be used, at the 
price of considerable computational complexity, to improve our knowledge of 
the in focus image. The idea, of course, is to use some generalized 
Gerchberg 3 algorithm to Iterate back and forth between the in-focus image and 
the out-of-focus diffraction pattern using 
o Known Fraunhofer diffraction laws, 
o Measured data in both domains, and 

o Imposed constraints in both domains (e.g. nonnegativity in the image 
plane). 

Such techniques are exceedingly powerful 1 * and would, no doubt, be useful here 
as well. 
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V. DEEP OBJECTS 


The hardest problems for any automatic, semiautomatic, or even human 
analysis Is the deep object. Such an object is deeper than the depth of focus 
of its hologram. It is never all in focus in any plane. Of such objects we 
ask the usual questions (shape, size, 3D orientation, 3D location), but 
finding the answers is quite difficult. 

The first problem which must be solved is that of separating in~f ocus 
images from out of focus artifacts. In any depth plane there are probably 
some of both. We suggest here an automated approach not as a final solution 
but as a starting place for more sophisticated analyses. 

We might begin by projection of the 3D scene into 2D. To do this we sort 
all x-y pixels by some focus criterion. That is at each x^,y^ we examine 
each discrete depth. We might estimate focus by brightness. Let the 
intensity at x^yj in the k th / depth slice be Iijk. We define 


P 


ijk 


max 



ijk 

Ij(k-D 


That is, Pijk, is the largest of the seen so far. By the time we 

have sorted through all N depth slices, PijN is a "projection" of the 3D 
image into 2D (i,j). All in focus pixels regardless of their depth are 
collected in one plane. 

The next step would be to reproject PijN back into 3D. The 3D 
reprojection of P^jjj is 


R ijk '{o 


ijk if *ijk 
otherwise 


ijN 



Thus Rijk should have no out of focus parts. The subsequent analysis must 
group non-zero R^j^ components Into likely particles, characterize them, 
etc. 


VI . PROSPECTS 

The sole objective of this paper Is to point out numerous opportunities 
to explore automated computer analysis of 3D particle fields obtained by 
holography. For shallow spheres the analysis is easy and only implementation 
Is needed. For other shallow objects some analogous work can be borrowed from 
the spherical case but much new work Is required. For deep objects, no work 
has been done. The critical observation, however, is that full automation 
appears to be within our grasp. 
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